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Abstract 

These are the proceedings of the workshop ’’HadAtomOl”, held at the Institut fiir 
Theoretische Physik, Universitat Bern, October 11 - 12, 2001. The main topics discussed 
at the workshop were the physics of hadronic atoms and in this context recent results in 
experiment and theory. Included here are the list of participants, the scientific program 
and a short contribution from each speaker. 
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1 Introduction 


The workshop “HadAtomOl” took place at the Institut fiir Theoretische Physik, Universitat 
Bern on October 11-12, 2001. It was the third in a series of workshops on bound states, in 
particular hadronic atoms: the previous ones were held in Dubna, Russia (May 1998) [1], and 
in Bern (October 1999) [2], The meeting was attended by about 50 physicists, and results were 
presented by 25 participants. 

The topics of the workshop included 

• Hadronic atoms, in particular their • Experiments 

Production DIRAC at CERN 

Interaction with matter DEAR at DAFNE 

Energy levels PSI (Pionic Hydrogen Collaboration) 

Decays Others 

• Meson-meson and meson-baryon scattering • decays 

The talks were devoted to recent experimental and theoretical progress in the investigations 
of hadronic atoms. Among the highlights of the workshop were the presentation of preliminary 
results of the DIRAC collaboration on the measurement of the lifetime of pionium, as well as 
the hrst measurement of kaonic nitrogen by the DEAR collaboration. 

The speakers have provided a one-page abstract of the presentations as well as a list of the 
most relevant references - these contributions are collected below. We also display a list of the 
participants with their e-mail, and the scientihc program of the workshop. 
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DIRAC (PS-212) experiment at CERN 

C. Santamarinal^ Q (on behalf of DIRAC collaboration) 

^ Departamento de Fisica de Particulas, Universidade de Santiago de Compostela, Campus 
Universitario Sur, 15782- Santiago de Compostela, Spain. 

The theoretical results in the low energy regime of hadronic physics have reached an accuracy 
level which should be tested with new and precise experimental measurements. In particular, a 
prediction of 2.9 ± 0.1 x 10“^® seconds [1] has been made for the lifetime of the ground state of 
pionium, the 7r+7r“ electromagnetic atom, within the framework of Chiral Perturbation Theory. 

DIRAC-PS 212 [2] is a hxed target experiment, currently running at CERN PS, which aims 
to measure the pionium lifetime with a 10% accuracy analyzing the low relative momentum 
region of 7r’''7r“ pairs spectrum {Q < 2bMeV/c) resulting from proton target collisions [3]. 

The DIRAC spectrometer is a double arm telescope with a central magnet. It is optimized 
to detect events with two oppositely charged particles with low relative momentum. The 
detectors provide an efficient particle identification to exclude events with proton, electron or 
muon contamination. 

The evidence of pairs from atomic break-up was first found in Summer 2000. The analysis 
of 2000 and 2001 data has given us around 5600 atoms mainly with a Titanium (1800 atoms) 
175 /rm and a Nickel (3600 atoms) 94 /rm targets. 

With the analyzed data a preliminary result, considering statistical errors only, has been 
obtained both with Nickel and Titanium data. The combined analysis of both samples gives 
us a lifetime of O.Olg;® x 10“^^s, however a detailed analysis of systematic errors is yet to be 
performed. 


Sample 
Target/Year 

Lifetime 

(xlO-i^s) 

Statistical 

error 

Ni 2000 


±34% 

Ti 2000/1 

5.471:^ 

±26% 

Ti+Ni 

'J-O-0.7 

±22% 


Table 1: Preliminary results with 2000 and 2001 data. Only statistical errors have been ana¬ 
lyzed. 


[1] J. Gasser, V.E. Lyubovitskij, A. Rusetsky and A. Gall, Phys. Rev. D 64 (2001) 016008. 

[2] B. Adeva et al, CERN-SPSLC-95-1 ; SPSLC-P-284 (1994). 

[3] L.L. Nemenov, Yad. Eiz. 41 (1985) 980; Sov. J. Nucl. Phys. 41 (1985) 629. 


"^Now at Basel University. 
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Spectrum of pious from tt+tt -atom breakup 
(Born and Glauber approximations) 

L. Afanasyev^, A. Tarasov^’^ and O.Voskresenskaya^ 

^ Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia 
^ Institute for Theoretical Phisics, University of Heidelberg, Philosophenweg 19, D-69120, 

Heidelberg, Germany 

Momentum and angular distributions in 7r“''7r“-pairs from the 7r’''7r“-atom (A 2 ,r) breakup 
at the electromagnetic interactions with target atoms have been considered in the Born and 
Glauber approximations. Exact analytical expressions for the spectra for the initial states with 
orbital quantum number / = 0 (nS states) have obtained in the hrst Born approximation. 




In the hgure the distributions over the relative momentum of pions p and angle 6 with respect 
to the total momentum of atom are shown for the initial states with n = 1 and n = 10 for the 
nickel target. 

The formulas for the discussed spectra for the initial ground states of A 27 r which takes into 
account muli-photon exchange processes have been obtained in the following approach. It is 
shown that due to the small size of A 27 r and its big mass the interaction with an unscreened 
target nuclear dominantly contributes to the multi-photon exchange processes. The amplitude 
of the atom interaction with unscreened Coulomb potential can be expressed in a closed an¬ 
alytical form in the Glauber approximation. Thus one can consider the target atom as the 
screened Coulomb potential for accounting the single-photon exchange processes in the Born 
approximation and as unscreened one for multi-photon processes. 
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Possibility to measure the energy splitting between 2S 

and 2P states in 7r+7r“ atoms 


L.L. Nemenov 

CERN, CH-1211, Geneva 23, Switzerland and Dubna (JINR), Russia 

The difference AEn between the energy of atomic nS and nP states of A 2 n includes the 
contributions of the vacuum polarisation AE^^, and of the strong interaction AEf^ effects, 
where AE^ ~ 2ao + 02 , and Oq and 02 are the S-wave tttt scattering lengths with isotopic-spin 
quantum numbers 0,2 [1,2,3]. 

The value of AE^^ is well known from QED calculations [1,2,4]. For this reason the 
measurement of AEn will give, in a model independent way, the value of 2ao -l- 02 [5,6]. 

The method for measuring AEn is presented [7,8]. The same method allows also to measure 
the energy splitting {AEn) for atoms consisting of tt and K mesons, and to obtain the value of 
2 ai /2 + 03 / 2 , where ai /2 and 03/2 are the S-wave nK scattering lengths with isotopic-spin 1/2 
and 3/2, respectively. 


[1] A. Gashi et ah, Nucl. Phys. A 628 (1999) 101. 

[2] A. Rusetsky, private communication. 

[3] G. Efimov et ah, Sov. J. Nucl. Phys. 44 (1986) 296. 

[4] A. Karimkhodzhaev and R.N. Faustov, Sov. J. Nucl. Phys. 29 (1979) 232. 

[5] G. Colangelo, J. Gasser and H. Leutwyler, Phys. Lett. B 488 (2000) 261. 

[6] M. Knecht et ah, Nucl. Phys. B 457 (1995) 513. 

[7] L.L. Nemenov, Sov. J. Nucl. Phys.41 (1985) 629. 

[8] L.L. Nemenov and V.D. Ovsiannikov, Phys. Lett. B 514 (2001) 247. 



Hadronic atoms: review of theoretical aspects 

J. Gasser, I. Mgeladze, A. Rusetsky and I. Scimemi 

Institute for Theoretical Physics, University of Bern, Sidlerstrasse 5, CH-3012, Bern, 

Switzerland 

In order to fully exploit the high-precision experimental data on hadronic atoms which 
is available at present and which will be provided in the future, it is imperative to design a 
theoretical framework that describes this type of bound systems to the accuracy that matches 
the experimental precision. The Deser-type relations which are then used to extract the “purely 
strong” scattering lengths from the measured values of the strong energy shift in the ground 
state, and the partial decay width into hadronic channels, contain isospin-breaking corrections 
which depend on the hne structure constant a and on the quark mass difference rud — mu- The 
following questions arise. 

i) The scattering lengths are calculated in QCD at uiu = rud- How is this theory related to 

a framework that includes isospin violating effects, SU{3)c x [Note that, in view 

of the announced accuracy 0.3 % in the future measurement of the energy shift of the 7r“p 
atom by Pionic Hydrogen collaboration at PSI, the question of the precise dehnition of 
the pure QCD limit is not an academic one.] 

ii) With a given dehnition of the pure QCD limit, how does one calculate the isospin-breaking 
corrections to the bound-state observables? 

The answers are as follows. 

ad i) In two-havour QCD, the quark mass niu = nid and the strong coupling constant may be 
hxed from the requirement that the pion and nucleon masses are equal - by convention 
- to the experimental values of the charged pion and the proton masses, respectively. In 
the effective theory of SU{3)c x U{l)em, the values of the strong LECs stay put [1]. 

ad ii) We have developed a general theory of hadronic atoms, based on a merger of non- 
relativistic effective Lagrangian techniques and ChPT [3-6]. This theory has been suc¬ 
cessfully applied to the calculation of the isospin-breaking corrections in the 7r'''7r“ atom 
decay width [3,4,6], and the energy shift of the 7i~p atom [5]. For related work, see [7]. 
Applications to other systems are in preparation [8]. 

[1] J. Gasser, I. Mgeladze, A. Rusetsky and I. Scimemi, in preparation. This reference reviews 
earlier work on the subject. We have profited from unpublished notes by H. Leutwyler, whom 
we thank in addition for informative discussions on the subject. 

[3] A. Gall, J. Gasser, V.E. Lyubovitskij and A. Rusetsky, Phys. Lett. B 462 (1999) 335. 

[4] J. Gasser, V.E. Lyubovitskij and A. Rusetsky, Phys. Lett. B 471 (1999) 244 (1999). 

[5] V.E. Lyubovitskij and A. Rusetsky, Phys. Lett. B 494 (2000) 9. 

[6] J. Gasser, V.E. Lyubovitskij, A. Rusetsky and A. Gall, Phys. Rev. D 64 (2001) 016008. 

[7] X. Kong and P. Ravndal, Phys. Rev. D 61 (2000) 077506 ; D. Eiras and J. Soto, Phys. Lett. B 
491 (2000) 101, and references cited therein. 

[8] M. Schmid, J. Schweizer and P. Zemp, work in progress. 
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Spectrum and decays of tt+tt and tt p atoms [1] 

V. E. Lyubovitskij 

Institut fur Theoretische Physik, Universitdt Tubingen, 

Auf der Morgenstelle If, D-72076 Tubingen, Germany 

In order to carry out the precision experimental tests of QCD we proposed and developed 
effective theory [2,3] for the description of measured characteristics of hadronic atoms. We 
have derived the expressions for the decay width r 2 ,r 0 of the 7r’''7r“ atom [2] and strong energy 
shift of the r~p atom [3] in the framework of QCD (including photons) by use of effective held 
theory techniques. The results obtained contains all terms at leading and next-to-leading order 
in the isospin breaking parameters a ~ 1/137 and — m^. On the basis of these formulae, 
a numerical analysis was carried out in Refs. [2,3] at order in chiral perturbation theory 
(ChPT). 

The 7r’''7r“ atom decay problem is completely solved [2]. Numerical analysis in ChPT based 
on recent update results for the S'-wave tttt scattering lengths Oq and 02 gives 

r = p-/o = (2.9±0.1) X 10-'^ s. (1) 

Numerical analysis of the Tr~p atom strong energy shift is in progress because the effective 
electromagnetic coupling constant /i is still unknown. 

The perturbative chiral quark model [4] was applied to estimate the electromagnetic O(p^) 
low-energy constants of the ChPT effective Lagrangian that dehne the electromagnetic mass 
shifts of nucleons and hrst-order (e^) radiative corrections to the ttN scattering amplitude [5]. 
It was found that the value of /2 constant is consistent with model-independent result obtained 
by Gasser and Leutwyler [6]. Using obtained value for constant /i the leading isospin-breaking 
correction to the strong energy shift of the 7r~p atom in the Is state was calculated. The 
result = —2.8 x 10“^ is comparable to a prediction based on a potential model for the ttN 
scattering [7]: = —2.1 x 10“^. 

Acknowledgements. This work was supported by the Deutsche Forschungsgemeinschaft 
(DFG, grant FA67/25-1). 


[1] Work done in collaboration with J. Gasser, A. Rusetsky, A. Gall [2,3] and Th. Gutsche, A. Fa- 
essler, R. Vinh Man [4,5]. 

[2] A. Gall, J. Gasser, V.E. Lyubovitskij and A. Rusetsky, Phys. Lett. B 462 (1999) 335; 
J. Gasser, V.E. Lyubovitskij and A. Rusetsky, Phys. Lett. B 471 (1999) 244; J. Gasser, 
V.E. Lyubovitskij, A. Rusetsky and A. Gall, Phys. Rev. D 64 (2001) 016008. 

[3] V.E. Lyubovitskij and A. Rusetsky, Phys. Lett. B 494 (2000) 9. 

[4] V.E. Lyubovitskij, T. Gutsche, A. Faessler and E.G. Drukarev, Phys. Rev. D 63 (2001) 054026. 

[5] V.E. Lyubovitskij, T. Gutsche, A. Faessler and R. Vinh Mau, Phys. Lett. B 520 (2001) 204. 

[6] J. Gasser and H. Leutwyler, Phys. Rep. 87 (1982) 77. 

[7] D. Sigg, et al, Nucl. Phys. A 609 (1996) 310. 
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Electromagnetic corrections to the lifetime of pionium 

G. Rasche 

Institute for Theoretical Physics, University of Zurich, Winterthurerstr. 190, CH-8057, 

Zurich, Switzerland 

We have made use of the Deser formula for the two-channel system as derived in [1], In¬ 
cluding a tiny correction for vacuum polarisation etc. in the numerical constant, we get from 
[ 2 ]: 

0.09474 

^ (fe) = 7-7^^ • (1) 

(ooc (fm))^ 

7r“ —>■ 7r°7r° transition element of the K-matrix of scattering 
It consists of a hadronic part aQ^(juo) and an electromagnetic 


Here aoc is the value of the tt" 
theory at the tt+tt” threshold, 
correction Aaoduo)- 


®0c — ®oc(ho) + ^floc(ho) • ( 2 ) 

Since the pion mass for the hadronic limit is very close to the neutral pion mass po if 
follows that the hadronic starting point for the true electromagnetic corrections is po- We have 
indicated this explicitly on the Ihs of (2). 

The value for is taken from [3]; they use aP — a? = 0.265(4). Taking into account 

that this value corresponds to a hadronic starting point with pion mass pc and is given in units 
of one gets (see [2]) 




4(ho) = -^(a (po) - a (ho)) = -0.170(3) fm . 


( 3 ) 


The electromagnetic correction Aaoc(ho) is calculated using a relativised two-channel Schro- 
dinger equation with two isospin invariant energy independent potentials V\r), / = 0, 2. They 
are constructed by htting an analytical ansatz to the two-loop hadronic ChPT scattering phase 
6\q), I = 0, 2 of [4], for both values po and Hc of the pion mass. Solving the coupled Schrodinger 
equation including the Coulomb interaction and the mass difference he — ho, as well as taking 
into account the results of [5], we arrive at 


Aaoc(ho) = 0.010(1) fm . 


Inserting these results into (1) we get 


This agrees with [3]. 


T = 2.92(15) fs . 


( 4 ) 


( 5 ) 


[1] G. Rasche and W.S. Woolcock, Nucl.Phys. A 381 (1982) 405. 

[2] A. Gashi, G.G. Oades, G. Rasche and W.S.Woolcock, |arXiv:hep-ph/0108116 , to appear in Nucl. 
Phys. A. 

[3] J. Gasser, V.E. Lyubovitskij, A. Rusetsky and A. Gall, Phys. Rev. D 64 (2001) 016008. 

[4] J. Bijnens, G. Golangelo, G. Ecker, J. Gasser and M.E. Saino, Phys.Lett. B 374 (1996) 210; 
Nucl. Phys. B 508 (1997) 263 and erratum ibid. B 517 (1998) 693; G. Golangelo, private 
communication. 

[5] M. Knecht and R.Urech, Nucl. Phys. B 519 (1998) 329. 
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Hadronic potentials from effective field theory 

E. Lipartia^’^’^, V. Lyubovitskij^’^ and A. Rusetsky^® 

^ Department of Theoretical Physics 2, Lund University, Solvegatan I 4 A, S22362 Lund, Sweden 
^ LIT, Joint Institute for Nuclear Research, 141980, Dubna, Russia 
^ HEPI, Tbilisi State University, University st. 9, 380086 Tbilisi, Georgia 

^ Institute of Theoretical Physics, University of Tubingen, Auf der Morgenstelle I 4 , D-72076, 

Tubingen, Germany 

^ Department of Physics, Tomsk State University, 634050 Tomsk, Russia 
® Institute for Theoretical Physics, University of Bern, Sidlerstrasse 5, GH-3012, Bern, Switzerland 

The problem of the evaluation of isospin symmetry breaking effects to the hadronic atom 
observables has been addressed in the framework of held-theoretical approaches, as well as in 
the potential models. At the same time, the results obtained within the potential model do not 
always agree with the results of held-theoretical approaches. The main reason of discrepancy 
between the predictions of both approaches is that, in general, the potential model does not 
take into account the full content of isospin-breaking effects in ChPT. Even if these predictions 
agree in a particular case (e.g., like the latest prediction of the potential model for the lifetime 
of 7r'''7r“ atom [1]), without an explicit matching to ChPT in the isospin-breaking phase the 
treatment of the problem in the potential model can not be regarded as systematic. In the view 
of the fact that the potential approach is widely used to evaluate the isospin-breaking effects in 
scattering process - to analyze the results of experimental measurements - it is very important 
to construct the potential which includes the full content of isospin-breaking effects. 

The constructive algorithm for the derivation of the hadronic potential from ChPT is con¬ 
sidered in detail in our paper [2]. Here we provide only a brief outline of the procedure. The 
derivation of the isospin-breaking part of the short-range strong potential is based on the uni¬ 
versality conjecture, which states that the relation between the bound-state characteristics and 
scattering amplitudes is the same in the held theory and in the potential model. The universal¬ 
ity conjecture provides us with the matching condition for the isospin-breaking part of hadronic 
potential, assuming the isospin-symmetric part of the potential is given. Further, the matching 
condition imposes rather loose constraints on the potential. This looseness rehects the freedom 
of the choice of regularization of pointlike singularities which emerges if the potential is derived 
from a local held theory. For this reason, the shape of the potential does not bear a physical 
information, and using this freedom, one may absorb all information about the isospin-breaking 
ehects in the couplings of the short-range potential which, in its turn, are determined from the 
matching condition [2]. 

Despite the unrestricted freedom in the choice of the shape of the short-range part of 
the potential, the full potential that we construct, reproduces, by construction, the threshold 
amplitude and the bound-state characteristics of hadronic atoms at the hrst non-leading order 
in isospin-breaking. 


[ 1 ] 

[ 2 ] 


G. Rasche, these proceedings. 

E. Lipartia, V.E. Lyubovitskij and A. Rusetsky, arXiv:hep-ph/0110186 
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Distribution of pious from breakup of pionium 

T.A. Heim ^. K. Hencken^, M. Schumann^, D. Trautmann^ and G. Baur^ 

^ Institut fiir Physik, Universitdt Basel, Klingelbergstrasse 82, CH-4056, Basel, Switzerland 
^ Institut fiir Kernphysik, Forsehungszentrum Jiilich, 52425 Jiilieh, Germany 

The experiment DIRAC requires in its analysis a reliable calculation for the number of pionium 
atoms produced in the target. At small relative momentum Q < 3MeV/c the observed momen¬ 
tum distribution shows an excess of tt+tt" pairs (attributable to ‘atomic’ pairs) as compared to 
the momentum distribution obtained from a fit at larger relative momentum where only ‘free’ 
pairs contribute [1]. However, extracting the number of pionium atoms in this way requires 
accurate (momentum) distributions of pions from the breakup of pionium as an input. 

In the framework of the semiclassical approximation [2], the distribution of the pions from 
the breakup, resolved with respect to magnitude and angular direction of the momentum of 
the out-coming pions, is given by 

where afi{b) denotes the amplitude for the transition from the bound initial state i to the 
hnal continuum state with (internal) relative momentum p of the pionium. Here, p = niredl = 
|(p_i_ — P-), with p±. denoting the momentum of the two pions in their center-of-mass frame. 
With the method outlined in [2] we can calculate these transitions into the continuum directly 
for arbitrary initial states. The azimuthal symmetry with respect to the beam axis reduces the 
angular distribution to a mere 6'-dependence where 6 denotes the angle between the beam and 
the out-coming tt’*' (or 7 i~, since the distribution is symmetric under the reflection 6 ^ tt — 6 in 
hrst order of the interaction, i.e., considering only breakup of the pionium due to one-photon 
exchange and taking into account the nearly perfect ; 2 -parity conservation [3,4], which becomes 
exact in the sudden limit). 

We have calculated and analyzed the differential cross sections for breakup of pionium on a 
Nickel target at projectile energy 5 GeV, for initial states n^s and n^p with n* < 4. We found 
that when integrated over the momentum, the angular distribution shows the two pions to 
be preferentially emitted perpendicular to the beam, as expected. This strongly non-isotropic 
distribution is to be contrasted with the isotropic distribution of the ‘free’ pairs not stemming 
from the breakup of pionium. 

When integrated over the angles, the momentum distribution for breakup from an initial 
state with principal quantum number Uj has its peak at p ~ (l/2nj) MeV/c, with a width of 
approximately Ap ^ (l/2nj) MeV/c. The peak amplitude scales as n^, as expected. Similar 
calculations for other initial states (d- and f-states reached by multiple scattering) and other 
target materials will be carried out. These momentum distributions will then be used in a more 
accurate extraction of the number of pionium atoms produced in the experiment. 

[1] L.G. Afanasyev et al, Phys. Lett. B 338 (1994) 478. 

[2] Z. Halabuka et al., Nucl. Phys. B 554 (1999) 86. 

[3] T.A. Heim, K. Hencken, D. Trautmann and G. Baur, J. Phys. B 33 (2000) 3583. 

[4] T.A. Heim, K. Hencken, D. Trautmann and G. Baur, J. Phys. B 34 (2001) 3763. 
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^ Institut fur Physik, Universitdt Basel, Klingelbergstr. 82, CH-fOSG, Basel, Switzerland 
^ Institut fur Kernphysik, Forschungszentrum Jiilich, 52425 Jiilich, Germany 


The total and excitation cross sections for the interaction of pionium with the target atoms 
are an essential theoretical input for the experiment DIRAC [1]. They are well known in the 
Born approximation and can be calculated to a very high accuracy including various corrections 
[2-4], However, Afanasyev et ah [5] showed that the total cross section in Born approximation 
overestimates the value in Glauber approximation by up to 14% for Z = 73 (Ta). This is 
a major setback for the prospect of calculating the cross sections to 1% as is required for the 
success of the experiment. To calculate the excitation cross sections we start with the transition 
amplitude between an initial state i and a final state / given by 

= / d^r'if}{r) [1 - exp iix{b,r))] ^pi{r} , 

where are non-relativistic hydrogen-like wave functions and x{b,r) = —{e/hv) /[<h(r’'/2) 

— ^{—r'/2)]dz with the screened potential $(r') = [Ze/r') X) exp(—a^r'). Unlike in the case 
of the total cross section, where the closure approximation leads to signihcant simplihcations, 
the excitation cross sections cannot be reduced to a one-dimensional integral in general. Instead 
a full numerical integration is required. Since the cross sections in Born approximation are well 
known, we only need to calculate the higher order corrections. Therefore we evaluate the 
integral 


Afi = 2% 


bdb[\afi{b)\ - |a^i(fe)| 


which converges much faster than the full cross section since the higher order corrections are 
most important for small impact parameters b. Figure 1 shows the results for a few selected 
transitions. One can see that even for a Nickel target {Z = 28) the higher order corrections are 
between 3% and 5%, while for large Z they reach up to 30%. Clearly, we need to take these 
results into consideration to reach the desired accuracy of 1% for the cross sections. 


Figure 1: Relative difference between 
Glauber and Born approximation for 
the excitation cross sections of ls-2p, 
2p-3d, and 2s-3p transitions as a func¬ 
tion of the target charge number Z. 


Z 



[1] B. Adeva et ah, CERN/SPSLC 95-1, SPSLC/P 284, Geneva (1995). 

[2] Z. Halabuka, et ah, Nucl. Phys. B 554 (1999) 86. 

[3] T.A. Heim, K. Hencken, D. Trautmann, and G. Baur, J. Phys. B 33 (2000) 3583. 

[4] T.A. Heim, K. Hencken, D. Trautmann, and G. Baur, J. Phys. B 34 (2001) 3763. 

[5] L. Afanasyev, A. Tarasov, and O. Voskresenskaya, J. Phys. G 25 (1999) B7. 
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Calculation of electromagnetic breakup of pionium 

T. Heim^, K. Hencken ^. M. Schumann^, D. Trautmann^ and G. Baur^ 

^ Institut fiir Physik, Universitdt Basel, Klingelbergstr. 82, CH-4056, Basel, Switzerland 
^ Institut fiir Kernphysik, Forschungszentrum Jiilich, D-52425 Jiilich, Germany 

The calculation of the electromagnetic excitation and breakup of Pionium in the target is 
an important input for the analysis of the experiment DIRAC at CERN. The cross section 
for excitation, as well as, breakup was studied within the semiclassical approximation for the 
dominant Coulomb part of the interaction (“scalar interaction”) in [1], In addition to coherent 
excitations, where the atom remains in the ground state, incoherent excitations have to be 
considered also. Summing over all excited states using the closure approximation one obtains 
a compact expression in terms of the inelastic scattering function Simi [2]. Here only the 
dominant Coulomb interaction is considered. The Simi were calculated within a Dirac-Hartree- 
Fock-Slater calculation. Simpler models like “antiscreening” and the “no-correlation limit” 
were found to be not adequate for the accuracy needed. The contribution of transverse photons 
was estimated to be at most 0.4% by making use of the long-wavelength approximation. 

Calculation of higher order effects are studied within the Glauber approximation and are 
discussed in [3]. This leaves us now with the calculation of the “magnetic interaction”, that is, 
the current interaction as the estimate show it to contribute almost 1%. In addition we want to 
study possible effects due to the relativistic treatment of the interaction between pionium and 
target [4]. The nonrelativistic Schrodinger wave function for the Pionium states could easily be 
replaced by the solution of the Klein-Gordon equation, but this gives only corrections of order 
Starting from the Klein-Gordon equation, we separate terms coming from the external field 
Ax, treating them as a perturbation. Using the Feshbach-Villars reduction to separate “large” 
and “small components”, we hnd the interaction between pionium and target to be 

ic r ^ ^ ^ ^ \ 

ieAxo + 7^— doAxo + VAx + AxV — -— (Axo ~ ^x) - ^cAxo- 

Besides the usual interaction terms, two new contributions are found; The interaction quadrat- 
ically in A is given by 

AU A _ a2 T2 _ ^ a.2 

^X^Xfj, — ^XO — ^ 2 ^ 0 ) 

which is suppressed by l/y^. The nonrelativistic treatment only leads to instead, which is 
not necessarily small. The second term is a Seagull diagram of the form <hcAxo with <I>c the 
Coulomb potential between 7r’''7r“. The same terms are also found starting from a diagrammatic 
approach, taking into account that Pionium consists of two equally heavy constituents. Explicit 
calculations of both the magnetic term, as well as, the Seagull diagram are performed. The 
Seagull diagram is found to contribute only to order 10“^; the magnetic terms are found to be 
less than 1%, in agreement with the estimates. 

[1] Z. Halabuka, et ah, Nucl. Phys. B 554 (1999) 86. 

[2] T. A. Heim, K. Hencken, D. Trautmann, and G. Baur, J. Phys. B 33 (2000) 3583. 

[3] M. Schumann et ah, contribution to this miniproceedings. 

[4] T. A. Heim, K. Hencken, D. Trautmann, and G. Baur, J. Phys. B 34 , (2001) 3763. 
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tiK atom: observation and lifetime measurement with 

DIRAC 


A. Lanaro 

(on behalf of the DIRAC Collaboration) 

CERN, CH-1211, Ceneva 23, Switzerland 

The detection of ttK atoms and the measurement of their lifetime provide a model in¬ 

dependent way to determine |ai/ 2 — 03 / 2 ], the difference between 1=1/2,3/2 S-wave nK scattering 
lengths. The nK scattering probes the SU{3)l x SU{3)r predictions of Chiral Perturbation 
Theory (ChPT) with u, d, s, and unequal mass kinematics, whereas the investigation of nn 
scattering checks the symmetry breaking part of SU{2)l x SU{2)ji with u and d quarks only. 

The nK scattering amplitude has been evaluated at next-to-leading order in standard ChPT 
[1]. On the experimental side, S'-wave nK phase shifts have been measured from threshold up 
to 1.3 GeV, using the scattering of moderately energetic charged kaon beams (2 < Ek < 13 
GeV). The results for the 1=1/2 and 1=3/2 S-wave nK scattering lengths are spread over a 
wide range: 0.13 < ai /2 < 0.24, —0.13 < 03/2 < —0.05, in units of inverse n mass. 

To determine the nK threshold parameters, the DIRAG Gollaboration proposes the same 
experimental approach as in the nn case, which consists in the detection of nK atoms [2]. 

A.„k are detected because they have a sizeable probability to break-up in traversing the 
target material where they are produced. The ratio of experimentally detected atomic 

pairs to the full number of produced atoms, yields a measurement of the break-up probability 
in a given target, and, thus, of the Attx lifetime. The relation between the ground state lifetime 
and the difference |ai /2 — 03 / 2 ! is given, in first order, by the Deser formula [3]. Using the GhPT 
prediction at leading order in Isospin breaking for ai /2 and 03 / 2 ] [1], the lifetime of the A^^-k 
ground state is predicted to be 4.8 ■ 10“^^ s. 

The Gollaboration envisages to use the present detection apparatus, as well as incoming 
particle beam and energy (24 GeV/c protons from the GERN PS). Only some new implemen¬ 
tations to the particle identification system would have to be considered to allow the simultane¬ 
ous detection of charged pions and kaons. Additional threshold Gherenkov counters of different 
type have been proposed [2]. Two identical detectors will employ a heavy gas radiator (Freon 
114), suitable to discriminate from K^ in the kinematic range of accepted A^^k- Another 
Gherenkov detector will employ Aerogel as radiator to discriminate K~^ from protons. 

In the summer of 2001 a test has been performed using the present Gherenkov detectors 
filled with SFe gas. The discrimination of pions from kaons was investigated by detecting the 
0(1020) — K^K~ resonance decay mode. During 2002, the Gollaboration will also perform 
tests of Aerogel samples with very low index of refraction (n=1.008), manufactured by the 
Boreskov Institute of Gatalysis at Novosibirsk. 

In three years of running time, DIRAG expects to detect a few thousands nK atomic pairs, 
which should lead to a determination of |ai /2 — 03 / 2 ] with 10-i-15% accuracy, and to a substantial 
improvement on the precision of the nn S-wave scattering length parameter |ao — 02 / 

[1] V. Bernard, N. Kaiser and U. Meissner, Nucl. Phys. B 357 (1991) 129. 

[2] B. Adeva et ah, the DIRAC Collaboration, CERN/SPSC 2000-032, 17 Aug 2000. 

[3] S.M. Bilenky et ah, Sov. J. Nucl. Phys. 10 (1969) 469. 
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Can one extract ms/m from ttK atom properties? 

H. Sazdjian 

Groupe de Physique Theorique, Institut de Physique Nueleaire, 

Universite Paris XI, F-91406 Orsay Cedex, France 
E-mail: sazdjian@ipno.in2p3.fr 

The lifetime and the 2p — 2s level energy splitting of the ttK atom are stndied. The lifetime 
is essentially dependent on the combination of the scattering lengths of the strong 

interaction process ttK ttK. This combination is independent, at leading order of chiral 
pertnrbation theory, of the qnark mass ratio mg/rh. The measurement of the ttK atom lifetime 
would therefore mainly provide information about higher-order (one-loop level) low energy 
constants. The energy splitting between the 2p — 2s levels is essentially dependent on the 
combination (2al^^ of the scattering lengths. This combination depends, at leading order 

of chiral perturbation theory, on the quark mass ratio mg/rh. This ratio is however multiplied 
by a term proportional to the deviation of the SU (2) x SU (2) quark condensate from the Gell- 
Mann-Oakes-Renner value. Therefore, if at the SU{2) x SU{2) level (tttt scattering amplitude) 
the chiral perturbation theory is close on experimetal grounds to the standard scheme of Gell- 
Mann, Oakes and Renner and of Gasser and Leutwyler, then the mg/m dependence of the 
2p — 2s energy splitting would hardly be detectable. On the other hand, if sizable deviations 
from the standard scheme are already present at the SU (2) x SU (2) level, then these are 
amplihed in the 2 p — 2s energy splitting by the quark mass ratio and experimental detection 
of the mg/m dependent term becomes possible. In any event, the measurement of the 2p — 2s 
energy splitting in the ttK atom would provide an experimental test for the standard scheme 
of chiral perturbation theory. In general, the measurements of the lifetime and of the energy 
splitting would allow the determination of each of the scattering lengths and 
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tzK scattering and dispersion relations 
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Chiral perturbation theory provides the low energy effective theory of the standard model 
that describes the interactions involving hadronic degrees of freedom. As a non-renormalizable 
theory, it requires additional low energy constants (LEG) that have to be introduced at each 
order of the loop- or momentum-expansion. In SU{3) at next-to-leading order, ten LECs 
Lj(p) enter the calculations. However, these quantities cannot be fixed by chiral symmetry 
only; experimental information is needed to pin them down. Using the Ki^ form factors and tttt 
sum rules provide accurate estimates for some of the LECs (e.g. Li, L 2 , L^) but leave the large 
Ac-suppressed LEC L 4 essentially undetermined. This is due to the fact that in such processes, 
contributions from L 4 accidentally are suppressed. Such a suppression does not occur in ttK 
scattering [1]. Therefore this process is suitable for a more reliable determination of L 4 and 
also provides estimates for Li, L 2 and L 3 . We use dispersion relations to fix these quantities [2], 
relying on analyticity, unitarity and crossing symmetry only, which are the suitable tools for a 
comparison of experimental data with ChPT. By combining fixed-t and hyperbolic dispersion 
relations for T^{s,t,u) and T~{s,t,u) and rewriting the chiral amplitudes, we find an appro¬ 
priate matching of the chiral and the dispersive representations of these amplitudes. Saturating 
the dispersive integrals with the available data [3], we find the following estimates 

Ll{mp) = 0.84 ±0.15, ^ 2 ^ = 1-36 ±0.13, 

Ll{mp) = -3.65 ±0.45, Ll{mp) = 0.22 ± 0.3, 

yielding an estimate for L 4 with improved reliability. 


[ 1 ] 

[ 2 ] 

[3] 


V. Bernard, N. Kaiser and U.-G. MeiBner, Phys. Rev. D 43, 2757 (1991); Nucl. Phys. B 357, 
129 (1991). 


B. Ananthanarayan and P. Biittiker, Eur. Phys. J. C 19, 517 (2001) |hep-ph/0012023 |; 
thanarayan, P. Biittiker and B. Moussallam, Eur. Phys. J. C (2001), to appear. 


B. Anan- 


P. Estabrooks et ah, Nucl. Phys. B 133, 490 (1978); D. Aston et ah, Nucl. Phys. B 296, 493 
(1988); D. Cohen et ah, Phys. Rev. D 22, 2595 (1988); A. Etkin et ah, Phys. Rev. D 25, 1786 
(1982) 
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Isospin violation in pion—kaon scattering 


B. Kubis and U.-G. Meifiner 

Forschungszentrum Julich, Institut fur Kernphysik (Theorie), D-52425 Jiilich, Germany 

Pion-kaon scattering near threshold is one of the cleanest processes to test our understanding 
of chiral dynamics in the presence of strange quarks. The existing determinations of the S-wave 
scattering lengths are, however, plagued by large uncertainties, such that high hope lies in the 
extraction of these parameters from ttK bound state experiments [1]. 

In order to relate the lifetime of ti~K^ (or 7 ;~^K~) atoms to one particular combination of 
the pion-kaon scattering lengths, one has to make use of a modihed Deser formula for the decay 
into the neutral tt^K^ channel which includes next-to-leading order effects in isospin breaking, 

rTrO^-O oc — CLq^ + e) (1 + -^) (1) 

(see [2] for an equivalent formula for the case of pionium). Here, e represents the isospin 
violating correction in the scattering amplitude 7 r~K~^ —> tt^K^ at threshold, while K is an 
additional contribution only calculable within the bound state formalism. 

We have evaluated the parameter e up to one-loop order, including leading effects both in 
the quark mass difference and the electromagnetic hne structure constant a = jF k [3] 

(see also [4]). While mu—rud effects are forbidden at leading order for tttt scattering and pionium 
and thus numerically negligible in comparison to electromagnetic effects, there are contributions 
linear in mu — md to the ttK scattering lengths which are of comparable size to the corrections of 
electromagnetic origin. At one-loop order, virtual photon loops induce a kinematical divergence 
at threshold (the Coulomb pole) oc l/|qin| which has to be subtracted in order to properly 
dehne the scattering length. Furthermore, photon loops contain an infrared divergence, which 
however vanishes at threshold. The corresponding radiative process tt~K^ —>■ which has 

to be included in order to dehne an infrared hnite total cross section yields a totally negligible 
contribution at threshold. 

The combined result for the scattering length, expressed in terms of relative corrections to 
the tree level isospin symmetric one, numerically reads 

ao(7r-A:+ ^ 7r°A:°) = ^ (ao^^ - 1 1 + 0^ + 0^ 

0{mu—md) 0{e^) 

+ (0.121 ± 0.009) + (0.008 ± 0.002) - (0.009 ± 0.008) 1 . (2) 

' -V- " ' -V- ' ' -V- ' J 

0(p4) 0{p^{mu-md)) 0{p^e^) 

The error ranges displayed for the next-to-leading order contributions are due to the uncer¬ 
tainties in the low-energy constants entering at that order. The value for the lifetime of the nK 
atom obtained using the “uncorrected” Deser formula, Tt^k = (4.8 ± 0.1) x 10“^® s, is slightly 
reduced to Tj^k = (4.7 ± 0.2) x 10“^® s. The additional correction K in (|ID however remains to 
be calculated. 


[ 1 ] 

[ 2 ] 

[3] 

[4] 


B. Adeva et al., CERN-SPSC-2000-032, CERN-SPSC-P-284-ADD-1, Aug 2000. 

A. Gall et al., Phys. Lett. B 462 (1999) 335; J. Gasser et al., Phys. Lett. B 471 (1999) 244. 

B. Kubis and U.-G. Meifiner, arXiv:hep-ph/010719E , Nucl. Phys. A, in print. 

A. Nehme and P. Talavera, arXiv:hep-ph/010729g . 
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In May 2001 the DEAR (DA<I>NE Exotic Atom Research) collaboration [1] performed the 
Erst measurement of an exotic atom (kaonic nitrogen) at the DA$NE collider of the Laboratori 
Nazionali di Frascati dell’lNFN. The measurement was performed with a cryogenic target (118 
K) filled with nitrogen at about 2 bar (p ~ 4.5p7VTp)- As detector, 8 CCD-22 were used. A 
total integrated luminosity of about 2 pb“^ was collected. After an optimization of the machine 
performances and of the DEAR degrader, an X-ray spectrum was obtained and analysed by 
performing a global £t. The global £t (with a y;^/NDOF = 0.95) gave a statistical signihcance 
of the two kaonic nitrogen lines measurable with DEAR of about 3 ct for the line at 4.58 keV 
(7 —>■ 6 transition) and of about 3.7a for the line at 7.6 keV (6 —5 transition), giving a global 
statistical signihcance of about 5 ct. These results compare favourable with the Monte Carlo 
simulation results with an yield of the two transitions of about An empty target-like 

measurement was as well performed. The background subtracted spectrum is presented in 
Figure 1. The experiment could therefore demonstrate the capability of creating and detecting 
exotic atoms using the kaons hux from DA<hNE. 

Next stage of the experiment is the measurement of the K-complex in kaonic hydrogen. 



Figure 1: The kaonic nitrogen background subtracted spectrum. X-ray lines are: 1 = A1 K^; 2 = Si 
K„; 3 = Sc K„; 4 = KN ( 4.58 keV); 5 = V K„; 6 = Cr K„; 7 = Cr Kp + Mn K^; 8 = Fe K„ + Mn 
Kp; 9 = Fe Kp; 10 = KN ( 7.62 keV); 11 = Cu K^; 12 = Zn Kq,; 13 = Zn Kp; 14 = Pb La 


[1] S. Bianco et ai, Rivista del Nuovo Cimento 22, No. 11 (1999) 1. 
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The quark condensate on the lattice 

P. Hasenfratz 

Institut fiir Theoretische Physik, Universitdt Bern, Sidlerstr. 5, 3012 Bern, Switzerland 


Chiral symmetry in QCD with Nf > 2 massless quark flavours is expected to be broken spon¬ 
taneously producing a non-zero {'ijj'ip) quark condensate generated through non-perturbative 
effects. A direct study of this condensate requires chiral symmetric lattice fermions (domain- 
wall, overlap, or the fixed-point formulations) which are expensive to simulate. The fermion 
condensate 


= — lim lim -— log Z 

mg^O V^oo dnir, 


( 1 ) 


is related to the low energy constant E of the tree level effective chiral lagrangean as {'ipip) = 
—Nf'E. A way to determine the low energy constant S is to measure the condensate in a small 
volume V in a fixed topological sector Q of QCD: Nf'Zrng,v,Q = —{'4’'4’)mgy,Q and use chiral 
perturbation theory to connect it to E defined at rUg = 0 and V = oo. Until now this is 
carried through in the quenched approximation only. Chiral perturbation theory predicts in 
this case [1] 


.. ^mg,V,Q - IQl/mqV _ E^u c 
rUg ~ 2\Q\ a2 ’ 


( 2 ) 


where \Q\/mgV is the contribution from the exact zero modes, while c/a^ (a is the lattice 
unit) is the ultraviolet divergent cut-off effect induced by the finite quark mass. Measuring 
Emq,y,Q in different volumes, eq. allows to And E. The bare low energy constant E is 
multiplicatively renormalized. The renormalization factor between E and Ej^(/i) is calculated 
by combining non-perturbative measurements and continuum perturbation theory. We obtained 
the result E;^(2GeV) = (270 ± lOMeV)^ [2], which, using the GMOR relation, corresponds to 
fnjig{2GeV) = (3.9 ± 0.4)MeV for the averaged u,d quark mass. This result is obtained at a 
resolution a ~ O.lSfm and no systematic continuum extrapolation is performed yet. The same 
applies also to the other measurements (they are consistent with our result) available at this 
moment [3]. One should keep in mind that these results refer to the quenched approximation 
{Nf = 0). It is unknown how far is the Nf = 0 low energy constant from that at Nf = 3 . The 
number above (which is close to the phenomenologically expected one) suggests that E(A^j) 
depends smoothly on Nf. This is not true, however for the condensate ('0'0) itself. The Nf = 0 
(quenched) QGD is not a healthy QFT and the tree level relation {'ip'ip) = —AtjE, unlike in full 
QGD, receives corrections which are even diverging in the chiral limit. As opposed to the low 
energy constant E, the condensate is not well defined in the quenched approximation [1]. 


[ 1 ] 

[ 2 ] 

[3] 


J.C. Osborn et ah, Nucl. Phys. B540 (1999) 317; P.M. Damgaard et ah, Nuch Phys. B547 
(1999) 305. 


P. Hasenfratz, S. Hauswirth, K. Holland, T. Jorg and F. Niedermayer, brXiv:hep-lat/0109007 . 

P. Hernandez et ah, JHEP 0107 (2001) 018; T. DeGrand, |arXiv:hep-lat /0107014 ; L. Giusti et 
ah, arXiv:hep-lat/0108007 . 
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In this talk I have reviewed the status of lattice calculations of the pion scattering lengths 
from a nonexpert point of view. The history of such calculations starts in the eighties/early 
nineties. I have split the discussion by reviewing hrst the calculations pre-2000 and then those 
post-2000. The most advanced calculation of last century is that of Fukugita et ah [1], where 
reference to earlier literature can also be found. This paper is also particularly relevant because, 
as far as I know, it is the only attempt at calculating the scattering length in the 1 = 0 channel 
on the lattice. Despite the huge effort and the remarkable results, various systematic effects 
were not under good control. 

A new generation of calculations by two Japanese collaborations have started in recent 
years [2,3]. Although the final results are not yet available, encouraging and very interesting 
preliminary results for the 1 = 2 scattering lengths have been presented at recent Lattice 
conferences. 

I have also discussed the recent results of a small Chinese collaboration that has exploited 
the use of improved actions to do the calculation on a very coarse and anisotropic lattice [4], 
and therefore perform it on a relatively small computer. I sincerely hope that this brave and 
encouraging calculation will challenge also smaller lattice groups to attempt such nontrivial but 
extremely interesting calculation. 

The agreement between experiments, CHPT, and lattice calculations for the tttt scattering 
lengths would represent a tremendous achievement for our understanding of strong interactions 
The CHPT calculation has reached a very high level of accuracy [5], and recent experimental 
results have confirmed the prediction [6]. Hopefully, the lattice calculations will also reach the 
same level of accuracy one day. 
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The quark condensate from i^e 4 decays 
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Roy equations [1] combined with chiral symmetry allow one to predict [2] the low-energy 
behaviour of the tttt scattering amplitude with high precision. The prediction is based on 
the hypothesis that the quark condensate is the leading order parameter of spontaneously 
broken chiral symmetry. This has been questioned in Refs. [3,4], where it is emphasized 
that experimental evidence for this scenario is not available. The recent result from the high 
statistics experiment E865 at Brookhaven [5] closes this gap, because those data allow 
one to determine the size of the leading term in the quark mass expansion of the pion mass. 
E865 conhrms [5] the hypothesis that underlies our prediction of the tttt S'-wave scattering 
lengths: more than 94% of the pion mass stems from the quark condensate. The generalized 
framework of SU{2) x SU{2) chiral perturbation theory developed in [3], that allows for a small 
or vanishing condensate, can therefore be dismissed. 

Note added'. In Ref. [6], the E865 data [5] have been analyzed without the use of chiral 
symmetry constraints. 
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Radiative corrections to semileptonic decays 

M. Knecht 
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Semileptonic decays of the light pseudoscalar mesons provide important informations on 
the non-perturbative aspects of the QCD vacuum (low-energy tttt scattering and low-energy 
constants Lj from decay), or on the CKM matrix elements Vud (pion /? decay) and 14s (-^€3 
decay). An accurate determination of these quantities requires high statistics experiments, as 
well as reliable theoretical descriptions. In particular, isospin breaking effects, arising both 
from rUu 7 ^ md and from radiative corrections, should be included. So far, this has not been 
done in a systematic way. 

A general framework allowing for a systematic treatment of isospin breaking effects within 
the chiral low-energy expansion has been introduced in [ 1 ], and applied to the simplest, 7^12 
and Ki 2 semileptonic modes. Recently, the channels have been treated as well [2] (for a 
preliminary discussion of the A '^4 modes, see [3]). 

In the presence of radiative corrections, the factorization between the hadronic and leptonic 
currents in semileptonic decays is lost. In the case of the Ki^ transitions, this entails that the 
two form factors f± not only depend on the momentum transfer between the two mesons, but 
also on a second, independent, kinematical variable. Furthermore, the necessity to include the 
decay mode with a soft photon emission, in order to obtain an infrared finite result, modifies 
the structure of the Dalitz plot density. Explicit expressions for the form factors at one 
loop with isospin breaking effects included have been obtained in [ 2 ], where existing expressions 
[4] for the Dalitz plot distribution have been corrected for several mistakes. 

The effects of isospin breaking on the extraction of I4s along the lines described in [5] have 
also been discussed and found to be small (see [ 2 ] for details). A reliable estimate of the two- 
loop strong interaction effects appears more than ever as being the crucial issue for a precise 
determination of I 4 s from decay [ 6 ]. 
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The framework used in Ref. [1] for studying K £4 decays is modified in order to account for 
isospin violating contributions, ^ rud and a ^ 0. The full kinematics, which requires three 
reference frames and hve variables introduced by Cabibbo and Maksymowicz Ref. [2] is changed. 
It is necessary to introduce a tensorial form factor, = ■j^Tpi^p 2 u, in the matrix element; 

However this tensorial form factor appears only in the 7 r’'“ 7 r“ channel, due to electromagnetic 
corrections. 

The decay rate is dTs = \ Vus P N{st,, si)J£,{st,, si,6t,,6i,(1)) x (is^(isz(i(cos6*^)d(cos6*z)(i0. We 

have a general formula for J 5 written in terms of the mass of the kaon which decays, the masses 
of the two pions in the hnal state and all the possible form factors, including the tensorial 
one. We explicitly calculate the form factors for at tree level using an effective Lagrangian, 
where not only the pseudoscalars but also the photon and the light leptons have been included 
as dynamical degrees of freedom Ref. [3]. 

The numerical values obtained for the decay rates(MeV) and phase space(MeV^) are shown in 
the table below. We compare to the values obtained by Bijnens and al. Ref.[4]. 


Decay 

(1) 7r"'"7r e~^Ve 

TT+TT 

(2) '!T°Tr°e~^h'e 

TT°Tr°fi'^iyu 

(3) 7r°7r e'^i'e 

7r°7r 

m±\t\ 

Tree level 

1297 

155 

683 

102 

561 

55 

Masses as in Rcf.[4J 
mu = md 

1292 

154 

681 

101 

561 

55 

e = 0 T = 0 
phase space MeV^ 

1.140 

0.146 

0.674 

0.099 

1.339 

0.197 

Physical masses 
mu / md 

1292 

153 

730 

no 

543 

51 

e / 0 T / 0 
phase space MeV^ 

1.140 

0.146 

0.674 

0.099 

1.307 

0.187 


In Ref.[4], the masses of the particles are (M^-, Mx) = (139.6, 493.6) MeV, (135, 493.6) MeV and 
(137, 497.7) MeV for the decays (1), (2) and (3), respectively. 

In case of no isospin violation, our Fortran routine reproduces the results of Bijnens and al. 
at better than 0.5 %. The interpretation of these numerical values leads to the result that the 
isospin violation, at tree level, is due to form factors for the channel and to phase space 
for the channel 7 r° 7 r“. The channel 7 r’'' 7 r“ is the only one which involves the tensorial form 
factor. At tree level, the effects of isospin violation and the existence of a tensorial form factor 
are numerically tiny in this channel. The calculations at one loop level are in progress. 
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At higher orders in the chiral expansion and when radiative corrections are included, a large 
number of unknown low-energy constants enters in ChPT. Lacking enough data, one often relies 
on estimates for these constants using some resonance Lagrangian. This approach was quite 
successful for the Lj at order [1], although one has to make sure that QCD short-distance 
constraints are met [2], Starting from the study of the low-energy and high-energy behaviors 
of the QCD three-point functions (VAP), {VVP) and (AAP), we have evaluated in [3] several 
0{p^) low-energy constants of the chiral Lagrangian within the framework of the lowest meson 
dominance (LMD) approximation to the large-iVc limit of QCD. In certain cases, values that 
differ substantially from estimates based on a resonance Lagrangian are obtained: 



00 

Cs 2 

CsT 

00 

Csg 

C 90 



Ay,P ‘2 

Ay,(p+g )2 


AA,(p+q)2 

LMD 

1.09 

-0.36 

0.40 

-0.52 

1.97 

0.0 

LAID 

-1.11 

-0.26 

-0.14 

0.040 








Set I 

-1.13 

0.0 

-0.096 

0.0 

Set I 

1.09 

-0.29 

0.47 

-0.16 

2.29 

0.33 






Set II 

1.49 

-0.39 

0.65 

-0.14 

3.22 

0.51 







In the table we have collected our estimates (LMD) for the low-energy constants Ci (even parity) 
from [4] and A* (odd parity) from [5], all in units of 10 “^/Fq^, and compared them to those 
obtained with a resonance Lagrangian, using different sets of parameters in the Lagrangian (set 
I and II). We have introduced the abbreviations Ay^p^ = A 2 — 4A3, = A 2 — 2 A 3 — 4 A 4 , 

Aa,p 2 = An + A 23 — 4A24, and AA,{p+q )2 = SAn -|- 4Ai7 -|- A23 — 2A24 — 4A25. 

The differences, in particular in Cgs, Cgo, Ay (p+q)2, and A^i (p+q)2, arise through the fact that 
QCD short-distance constraints are in general not correctly taken into account in the approaches 
using resonance Lagrangians. Whereas in the case of the vector form factor of the pion our 
result for a certain combination of low-energy constants at order p® agrees with the one [6] from 
the resonance Lagrangian, our estimate for one of the counterterm contributions in the form 
factor A in the decay vr —>■ ez/gq is by a factor of 5 smaller than the resonance estimate in [7]. 

Finally, in Ref. [3] we also briefly discussed how to systematically improve our approach. 
For instance, for the pion-photon-photon transition form factor it is necessary to include a 
second vector resonance (LMD-I-V) in order to correctly reproduce the CLEO data [8] for the 
form factor at high with one photon on-shell. In this way one obtains the new estimate 
Ay^-P~^^ = —1.36 X 10“"^/Fo, i.e., about 20% away from the LMD estimate given in the table. 
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The new pionic hydrogen experiment at PSI 
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A new high-precision determination of the strong-interaction shift (ei^) and broadening 
(Fis) of the ground state in pionic hydrogen (ttH) has been started at PSI [1]. This constitutes 
a direct measurement of the ttN scattering lengths aF and a“ and is an important test of the 
methods of chiral perturbation theory. 

The experimental approach is to determine the energies of x-rays emitted by radiative 
transitions to the Is-state in pionic hydrogen. This state is both shifted in energy because of 
the strong interaction between proton and pion and broadened because its lifetime is shortened 
due to the additional decay channels. The energy of the np-ls transitions in vrH is in the range 
of a few keV and should be determined with an accuracy in the 20 meV range. This can be 
achieved only with a crystal spectrometer. 

A previous experiment at PSI [2] measured the 3p-ls transition in pionic hydrogen and 
deuterium (ttD), combining a highly efficient stopping arrangement (cyclotron trap) for the 
pion beam and a cylindrically bent Bragg crystal with newly developed CCD (Charge Coupled 
Device) detectors. While an accuracy of 0.5% was reached for ei^, the error for Pi^ (8%) was 
an order of magnitude higher, which makes it hard to determine both a"*" and a“ from the ttH 
measurement alone. The ttD shift measurement can be used to improve the situation, but this 
necessitates a good theoretical understanding of the 3-body-problem in the deuterium system. 

The goal of this experiment is to reduce the uncertainty of Pi^ to ~1% and that of ei^ to 
0.2%. This can be achieved by using an improved cyclotron trap, spherically bent crystals and 
a new large-area CCD for a higher x-ray rate and better background suppression. 

During the hrst production run the is transition was measured at three different pressures 
(4 bar, 28 bar and liquid hydrogen) to investigate the possible influence of mechanisms that 
depend on the collision rate between vrH atoms and target molecules. One process is the 
formation of molecular states (tt^H -|- II 2 —[('7r“pp)p]ee) that could skew the ei^ measurement. 
The other is Coulomb deexcitation (7r“H(n) -|- H 2 —7r“H(n/) -|- H 2 ), a non-radiative process that 
leads to a Doppler broadening of the line. 

A preliminary analysis of the data (Table 1) does not show a pressure dependence of the shift 
eis. Furthermore, the averaged result does not deviate from the earlier result [2] which gave ei* 
= 7108 ±13 ±34 at 15 bar. The width Pi* (including the Coulomb deexcitation) appears to be 
pressure dependent. This will be investigated in the next stage of the experiment, starting in 
Spring 2002. 


Target pressure 

eis [meV] 

Pi. [meV] 

4 bar 

7082 ± 31 ± 15 

973 ± 75 ± 10 

28 bar 

7137 ±18 ±40 

969 ± 26 ± 10 

liquid (~700bar) 

7095 ± 25 ± 25 

1052 ±58 ±10 


Table 1: Preliminary fit results for the Summer 2001 data (value ±stat. ±sys.) 
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Atomic cascade in hadronic atoms 
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The atomic cascades in exotic hydrogen atoms x~p {x~ = , 7i~, K~, p) have been studied 

using a new Monte Carlo code that calculates the evolution of the x~p kinetic energy dis¬ 
tribution from the very beginning of the cascade. The cascade model is based on new cross 
sections for the collisional processes: Stark mixing, elastic scattering, nuclear absorption during 
collisions. Coulomb deexcitation, and external Auger effect. For the final part of the cascade 
(n < 5), Stark mixing, deceleration, and nuclear absorption during collisions for x~p scattering 
from hydrogen atoms have been computed in the close-coupling approximation [1-3]; this fully 
quantum mechanical framework takes threshold and strong absorption effects into account. For 
n > 8 , a classical-trajectory Monte Carlo model [3] for exotic atoms scattering from molecu¬ 
lar hydrogen has been used to calculate the cross sections for Stark mixing, elastic scattering, 
and Coulomb deexcitation. The semiclassical (eikonal) approximation is used for values of n 
between the quantum mechanical and classical domains. 

The predictions of the detailed kinetics cascade model have been confronted with the ex¬ 
perimental observables: X-ray yields, cascade times, and kinetic energy distributions [4]. The 
case of pionic hydrogen is of particular interest. The new PSI experiment [5], whose goal is to 
determine the strong interaction width F 15 of pionic hydrogen on the level of 1 %, relies crucially 
on a good understanding of the kinetic energy evolution during the atomic cascade because the 
measured nP IS' line profiles must be corrected for the Doppler broadening (see Fig. |l|). 



Figure 1: The calculated energy profile of the K /3 line at 
3 bar in pionic hydrogen with Doppler broadening taken 
into account (solid line) in comparison with the natural 
line shape assuming F 15 = 0.87 eV (dashed line). 
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Pion-nucleon analysis 
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There is information of the pion-nucleon interaction in different charge channels available 
at low energy from a variety of experiments: 1) cross section information, 2) analyzing powers, 
3) forward D'^ amplitude, 4) pionic hydrogen level shift and width measurements. Over the 
years such information has been analysed with different techniques, e.g. 

1. Karlsruhe fixed-f analysis KH78/80 [1], partial wave dispersion relation analysis KA.84 
[2] and partial wave relation analysis KA.85 [3]. 

2. VPI-GWU analysis with fixed-f constraints, accessible through the SAID facility 
http://gwdac.phys.gwu.edu/ 

3. Gashi et al. analysis at very low energies with an effective range expansion [4], 

4. Low-energy analysis exploiting forward dispersion relations [5]. 

A major problem in extracting hadronic quantities from observables is the question of the 
treatment of the electromagnetic effects. In low-energy scattering the formalism of Tromborg 
et al. [6] is often used, but there are effects not taken into account by that method. Also, 
for pionic hydrogen, determining the hadronic scattering lengths requires a careful analysis of 
electromagnetic corrections [7,8]. For the low-energy scattering region a 0{p^) calculation in 
chiral perturbation theory including all electromagnetic effects to that order has recently been 
published [9]. 

It turns out, especially in the analysis of sensitive quantities like the S-term, that in addition 
to the electromagnetic corrections also the high partial waves are important. Therefore, it is 
well motivated to start a new analysis in the spirit of Karlsruhe [10]. 
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In the Standard Model there are two sources of isospin violation. Strong isospin violation is 
driven by the quark mass difference (m^ — m^) whereas electromagnetic (em) isospin violation 
is due to the different quark charges. To extract the strong contribution, which leads to small 
effects of the size (m„ — mrf)/AQCD ~ 2%, one thus has to be able to simultaneously treat 
both sources. In the two flavor sector of QCD, pion-nucleon scattering is best suited to study 
isospin violation since the isovector operator {mu — md){uu — dd) contributes at the same 
order as the isoscalar chiral symmetry breaking one, {mu + md){uu + dd). In the isospin limit, 
the ttN scattering amplitude decomposes into isoscalar and isovector parts. This allows to 
formulate a variety of triangle relations between measurable channels that vanish if isospin 
were an exact symmetry, see e.g.[l]. As a first necessary step, one has to consider the isospin 
limit. In the framework of heavy baryon CHPT, the S- and P-wave phases can be described 
precisely in a fourth order calculation [2]. Fixing the low-energy constants (LECs) for pion lab 
momenta between 50 and 100 MeV by htting to the existing ttN partial wave analyses (KH, 
VPI/GWU, Zurich), one can predict threshold parameters and the phases at higher energies. 
Most importantly this calculation shows convergence of the chiral expansion. The inclusion of 
virtual photons (denoted as 7 *) in the effective ttN Lagrangian is straightforward, for a detailed 
discussion see e.g. [3]. In [4], a complete third order calculation was performed including one- 
and two-photon exchange diagrams, photon loop corrections, soft photon radiation and three em 
LECs, which parameterize the effects of hard photons and vacuum polarization. A £t to all data 
(with the exception of the ones generally considered inconsistent) below pion lab momenta of 
100 MeV allowed to determine the strong and em LECs and to uniquely separate the hadronic 
and the em amplitudes. A very pronounced difference in the hadronic amplitude for elastic 
scattering n~p n~p compared to the standard PWA’s was observed. This difference could be 
traced back to the inclusion of non-linear TnrNN'j* couplings as demanded by chiral symmetry. 
Such effects have to be included and thus a revision of the commonly used em corrections like 
e.g. the NORDITA approach is called for. In the S-wave triangle relation between elastic 
scattering, > n^p, and charge exchange, 7i~p —> n^n, one finds strong isospin violation of 

—0.75% in the low-energy region, an order of magnitude smaller than reported in the literature 
[5] but consistent we the expected size of isospin violation. It is now of utmost importance 
to map these result obtained in chiral perturbation theory on a dispersive representation to 
analyze also ttN data at higher energies. 
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